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We show that the ideal gas of Bogolyubov quasi-particles (IGBQ) of 2D Coulomb Bose gas 
provides c ~ T 4 and c ~ T (from plasmon and free particle parts of quasi-particles, respectively) 
temperature dependencies of the specific heat. This can in a framework of Drude model interpret 
the observed low-temperature heat conductivity k ~ T 3 ' 6 and k ~ T of pseudogap quasi-particles 
of cuprate superconductors. Vanishing of free quasi-particles of IGBQ in the transition to plasmons 
determines the downturn behavior of k and clarifies the origin of decoupling of charge carriers 
with heat. Near to pseudogap boundary Bose gas transforms into Fermi gas of holes or electrons. 
Obtained from mixture of Bose and Fermi gases Lorentz ratio values of Wiedemann - Franz law 
(WFL) are close to observed non-Fermi liquid ones. The observed temperature dependence of 
entropy can be result of the specific heat contributions from these both gases. 

PACS numbers: 74.72.-h, 74.20.Mn, 74.25. Fy, 74.25.Bt 



In the seminal paper [jj on investigation the low- 
temperature WFL for heat transport, Hill with col- 
laborators have experimentally studied the properties 
of normal state quasi-particles in the electron doped 
Pr2- x Ce x CuOi copper-oxide. It was shown by Hill et al. 
that they failed the Landau Fermi liquid theory (LFLT) . 
In the high external magnetic field the superconductiv- 
ity was suppressed and the heat conductivity of quasi- 
particles had a linear k ~ T law at a decreasing the tem- 
perature up to some value and then undergone k ~ T 3 
downturn behavior. The latter was called in the litera- 
ture a downturn of heat conductivity. Next important 
result of [l| was Lorentz ratio did not equate to Som- 
merfeld's value. These two observations were regarded 
as evidences of breakdown of LFLT in cuprate supercon- 
ductors. 

The downturn was observed in Ref. [2| and in 
the heavily over-doped non-superconducting compound 
La2~xSr x CuOi at low temperatures. The violation of 
WFL in the normal state of Bi2+ x Sr2- x C'uOe+s copper- 
oxide in the vicinity of metal-insulator transition was also 
detected in the recent paper 

Now it is widely believed [|[ that the phenomenological 
explanation of downturn is in the decoupling of phonon 
thermal bath with current producing charge carriers [4] 
and therefore, has no relation to fundamental properties 
of high-T c superconductors. In contrast, Hill et al. have 
supposed that the downturn may be intrinsic for copper- 
oxides. It would be interesting to find the experimental 
revealing of the downturn effect of heat conductivity in 
the another observable quantity, if it would be intrinsic, 
and understand the general underlying physics. 

It seems the experiment gives us this opportunity. 



Namely, no T dependence observation 0, H| of the elec- 
tric conductivity at low-temperatures allows to explicitly 
connect in the framework of Drude model the tempera- 
ture dependence of heat conductivity k with one of the 
specific heat c. 

The existing observed results on c of cuprate mate- 
rials are as follows. The normal state electronic spe- 
cific heat of high- T c superconductors Y Ba,2CuzO§j rX and 
La2- x Sr x Cu04 was experimentally investigated in [j| 
and [6(, respectively. Its magnetic field dependence of 
Yo,g,Cao.2Ba2Cu30e+ x compound has considered in 0]. 
The existence of superconductivity at low temperatures 
does not allow to extract the T dependence of the nor- 
mal state c. However, the normal state properties are 
clearly seen [5] in the entropy S. For instance, the dop- 
ing dependence of <S in the Y Ba2Cu^OQ+ x coper-oxide 
was displayed in It was also shown in [j| that for 
under-doped material the low-temperature dependence 
was S ~ T n with n > 1, while S ~ T was obtained for 
optimal doping compound. Non- linear T dependence of 
the entropy was related to isolating state Q of under- 
doped materials. 

One might assume that k ~ T 3 6 and S ~ T n have the 
same underlying physical origin and could be expressed 
through T dependence of the common specific heat c. In 
this case downturn of k would have a physical meaning 
and correspond to isolating state, which is consistent with 
the decoupling mechanism Q of thermal phonons with 
current carrying particles. 

In the present paper we follow Ref. 0] that high-T c 
superconductivity (HTCS) and pseudogap (PG) regime 
are described by mechanism of the single particle 2D 
Coulomb interacting bosons, in which the normal state 
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of superconductors is IGBQ. The physical origin of 
bosons is in the bosonization of 2D fermions below PG 
boundary due to coupling of their spins with statistical 
(anyon) magnetic field. According this mechanism in the 
deeply under-doped regime and low temperatures dom- 
inate bosons, while close to PG boundary there occurs 
the decoupling of spins and statistical magnetic field and 
dominate fermions. Therefore, the normal state is mix- 
ture of bosons and fermions. In 0] we have used the 
Bogolyubov approximation for 2D Coulomb Bose gas, 
which is valid only in the high density limit of gas. Being 
applied to gas from neutral atoms, this approximation 
has validity for weakly dilute densities. However, as it 
was shown in Ref. [9], it can give good results also for 
superfluid He II. Thus, we can assume that our results 
for bosons are applicable for intermediate densities of 2D 
Coulomb Bose gas. 

We do the calculation of the specific heat of IGBQ. 
For ideal Bose gas the thermodynamic free energy has 
the form [H: 



F = T 
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The energy scale is implied for temperature T. In the 
energy of quasi-particles [8j e(p) = (ap + (p 2 /2m) 2 ) 1 / 2 
ap is the square of plasmon energy with a = 2ire 2 nh/m, 
where e and n are the charge and 2D density, respec- 
tively, of true particles with mass m, and p is momentum 
of quasi-particles. 

The integral in Eq. (TTJ) over the absolute value of mo- 
mentum p we write in the form: = J 9 + . Here q 
provides at T = the maximum of the momentum dis- 
tribution function for true non-condensate particles [llT ] 
and separates the plasmon and kinetic energy contribu- 
tions to energy of quasi-particles. 

If for the first (second) integral we introduce the vari- 
able x = (apf/ 2 /T (x = p 2 /(2mT)) and then take into 
account that (aq) 1 / 2 = q 2 /(2m) we obtain 
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Here S is the area of 2D system. 

For the low temperatures, when (aq) 1 / 2 /T ^> 1, we 
can replace the upper limit of first integral in Eq. ([2]) 
by infinity sign. In this case the second integral becomes 
zero. For the high temperatures we have (aq) 1 / 2 /T 1 
then the lower limit of second integral can be putted zero 
and zero becomes the first integral. 

Following the scheme of calculation [l(|, we find the 
expressions for the specific heat divided to S 
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for ideal gas of plasmons and 
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for ideal gas of free quasi-particles. Here A = 
-(l/4)f™x i dx/(e x - 1) = -(l/4)r(5)C(5) and B = 
-J™xdx/(e x -l) = -r(2)C(2), where T(y) and ((y) are 
Gamma and the Riemann Zeta functions, respectively. 

The normal state heat conductivity k is described by 
WFL for free quasi-particles (charged bosons) . Its charge 
conductivity a, at low-temperatures, can approximately 
be obtained from the Drude model. The gas of plasmons 
does not care the charge therefore, is isolating. However, 
to determine the approximate downturn temperature Td 
we can formally extend the WFL to the region of this 
isolating gas. The transition from electrically conducting 
into non-conducting case will be described as transition 
from finite to zero (as it should be for isolators) charge 
conductivity. 

So, the estimate of downturn temperature Td is deter- 
mined from equation 



Kl(T d )/(7l = K 2 (T d )/cr 2 . 



(5) 



Here the Drude model gives the expressions of heat con- 
ductivities for gas of plasmons K\(Td) = (l/2)ci(Td)v\li 
and gas of free quasi-particles n 2 (Td) = (l/2)c2(Td)v2h- 
Corresponding mean velocities, free penetration lengths 
and charge conductivities are i>i.2, h.2 and cri^, respec- 
tively. 

The solution of equation Eq. ([5]) is 
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One can show that (ma 2 ) 1 / 3 = 2(2 / r 2 ) 2 / 3 Ry, where r s is 
the mean distance between particles in Bohr radius, as 
(as = H 2 /(me 2 )), units and Ry = me A /(2h 2 ) is Rydberg 
energy. 

We use the density of charg es n ab = 6.747- 10 14 • t ■ 
1 /cm 2 per 1 cm 2 area with doping fraction t of charge per 
Cu atom. It is almost the same for all cuprates, which is 
result of equivalent spacing constants [l2l | for elementary 
structural cells of C11O2 current carrying a — b planes. 
We rewrite this density through units then r s will be 
the function of t. Td has a form: 



T d « 0.114 



\G2Vlh 



t 2 ' 3 Ry. 



(7) 



Eq. ([7]) is an estimate of Td- We do the numerical 
calculation of the specific heat for arbitrary T to precisely 
obtain Td- 

The expression for the specific heat is 
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FIG. 1: The specific heat c/z (in q 2 /(8irh 2 ) units ) as function 
of z. 



where z = T/{aq) 1 / 2 and y = p/q . The result of nu- 
merical calculation of c/z (in q 2 /(8irh 2 ) units ) as func- 
tion of z is plotted in Fig. 1. This quantity equals to 
(Airh 2 /m)c/T. In the assumption that heat conductivity 
is described by Drude model, k = mv 2 ca / (2ne 2 \ and a 
does not depend from T (as was observed in 0, [|[), it 
displays temperature dependence of k/T (we note that 
for degenerate Fermi liquid and 2D Coulomb Bose gas v 
is function only of concentration n). Experiments [2|, H[ 
provide unique opportunity when only the specific heat 
c describes T dependence of k and we can compare this 
result for k with one for entropy S. 

From Fig. 1 we see that downturn of c/T oc- 
curs at (ag) 1 / 2 /2. Calculation gives (aq) 1 ^ 2 = 
(167re 2 nft/(3(2m) 1 /2))2/3 _ 8(2 1 /2/(3r 2 )) 2 / 3 i?y. Ex- 
pressing r s through t one obtains (aq) 1 ^ 2 — 0.737 -t 2 ^Ry 
and exact downturn temperature T^, Eq. [[?])■ with nu- 
merical factor 0.368. 

We estimate quantities at t ~ 0.1. From expres- 
sion for (ag) 1//2 one derives (aq) 1 ^ 2 jhs ~ 10 4 if, there- 
fore, the downturn temperature of the specific heat takes 
place at 10 4 -ftT. Observed downturn of heat conductiv- 
ity [E i] is T d}ex p/k B ~ 0.1K. Using Eq. © we find 
(T1V2I2 1 ~ 10~ 12 . One can assume v\l\ ~ 1^2, 
i.e. the product of mean velocities and free penetration 
lengths for plasmons and free quasi-particles has the same 
order of magnitude. Hence, we obtain cri/a% ~ 10 -12 . 
This result reflects the above supposed transition from 
charge conducting gas from free quasi-particles to iso- 
lating one from plasmons, i.e. metal-isolator transition. 
Therefore, the downturn of heat conductivity might be 
result of metal-isolator transition. It is worth here to note 
that the infinitesimal ratio of charge conductivities maps 
large scale of the specific heat temperatures of IGBQ into 
low ones of heat transport. We will discuss the metal- 
isolator crossover of cuprates within the single particle 
boson mechanism of HTCS elsewhere. 

We calculate the Lorentz ratio for free quasi-particles 



of IGBQ. The WFL, k/(ctT) = L , of heat transport for 
quasi-particles of LFLT is expressed through Lorentz ra- 
tio L by Sommerfeld's value L = (tt 2 /3)(k B /e) 2 . We 
note that this value of Lo corresponds to the three di- 
mensional (3D) case. A simple calculation shows that 
it does not change and for 2D. For the IGBQ we as- 
sume that the mean kinetic energy of quasi-particles is 
mv 2 /2 = {aq) 1 / 2 (because as was pointed out above at 
T = the most part of non-condensate particles has 
this energy). Then the WFL, n/{aT) = mv 2 c/(2ne 2 T), 
with c determined from Eq. Q for free quasi-particles 
of IGBQ has a form k/(ctT) = 3.106 • L /t 1/3 . This ex- 
pression corresponds to WFL of bosons when there is no 
mixture of bosons with Fermi quasi-particles of LFLT. 
However, at the increasing of concentrations of holes or 
electrons t to the direction of PG boundary there appear 
the fermions in the PG region Q and this mixture oc- 
curs. Therefore, at T — > we can phenomenologically 
write the expression n/{aT) — L with 



L = Ln 



3.106 / t\ t 

V-'TJ + Z 
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which takes into account the transformation of L from 
one of free quasi-particles of IGBQ to one of LFLT quasi- 
particles when concentration t tends to (but below) crit- 
ical concentration t c ~ 0.19 (l3l |. 

We do a comparison of our L/Lq with the experimen- 
tal one of Q. In Fig. 7 (a) of [|[ the data of L/Lq for 
several copper-oxides are plotted as function of t. This 
dependence of t is almost universal for all cuprates. We 
use results for Bi2+ x Sr2- x CuO§+s compound and ob- 
tain for t: 0.143, 0.17 and 0.185 the theoretical values 
of {L/L )theor- 2.222, 1.485 and 1.117, while the exper- 
imental ones of (L/Lo) exp are 3.1 ± 0.7, 1.2 ± 0.2 and 
1.05 ± 0.1, respectively. 

Now before discussion of our result on entropy, we 
would like to outline the some observed qualitative fea- 
tures [1, H, [l3[ of the HTCS specific heat and entropy. 
They are as follows. The increment coefficient 7 = c/T 
of the specific heat in the HTCS-to-normal-state phase 
transition (PT) point has weakly apparent and washed 
out peak in the deeply under-doped site. The peak of 
7 and its form become higher and sharper, respectively, 
when increases t. Peak is maximal (with the sharpest 
form) at t c and then does not almost change the form for 
t > t c . There is no influence of the external magnetic 
field on the normal state c. It is observed the loss 0] 
of the normal state entropy S in cuprate when doping 
changes from optimal doping to under-doping, which is a 
result 5] of changing of S ~ T into S ~ T™ with n > 1. 
And last observation is 7 does not depend from doping 
for t > t c . 

In general, the weakly apparent peak of 7 is attributed 
[i~c| to first order PT, which is close to second order one. 
Sharpening of 7 at t c might characterize the transition 
of PT into second order type, where conventional super- 
conductivity scenario with Cooper pairs becomes actual. 
This hypothesis was pointed out in 
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No dependence of the normal state c from the magnetic 
field might be result that quasi-particle energy e(p) of 
IGBQ is independent of magnetic field. To outline it we 
follow the argument of Abrikosov, given in the book [ijj], 
for Cooper pair quasi-particle energy ec P (p), which is 
independent of magnetic field. He has shown that ec P (p), 
which is result of canonical Bogolyubov transformation, 
is function of transformational scalar quantities u p and 
v p . The latter ones in the magnetic field with gauge 
V -A(r) = can be function of only zero scalar product k- 
Aj: = 0, where k and A^ are Fourier transforms of vector- 
coordinate r and vector-potential A(f), respectively. 

The calculation of the specific heat Cf of 2D gas of 
LFLT quasi-particles gives the same linear T dependence 
as in Eq. |4]), but with ratio cf/c2 = 2. This is due to two 
directions of spins for fermions instead of bosons. From 
Eq. ((4]) it is seen that cf is concentration independent. 
Therefore, observed no t dependence of 7 for over optimal 
doping of copper-oxides proves the 2D nature of cuprates. 

We find the expression of the normal state entropy 
S. ft is obtained through the calculation of integral 
S = \c/Ti)dTi . Our goal is to consider again the PG 
region, thus T < T*, where T* is PG boundary tempera- 
ture. Typical experimental scale of T* is T*/k B ~ 10 2 K, 
while (aq) 1 ! 2 jks ~ lO 4 ^, therefore, in Eq. (jHJl one can 
assume z < 1 and use Eq. §3§ as approximate expres- 
sion of the specific heat of IGBQ. For arbitrary t and 
T (in the intervals of t < t c and T < T*) we write the 
phcnomcnological expression 



c_ _ £1 / _ _5T 
T ~ T \ 4T* 



cfT_ 



(10) 



Except T/T* both terms of this formula should also 
been include the term 7r£ 2 (i)/(7!"£ 2 (£ c )) with the same 
sign and numerical factors as for T/T* . Here 7r£ 2 (i) is 
area Q of effective spin fluctuations of bosons at con- 
centration t and T = 0, where bosons reveal fermion 



property and their spins interact with spins of Anti- 
Ferromagnetic environment atoms of cuprates. From 
Fig. 4 of paper H we find that 7r£ 2 (t — > t c ) — > 00 and 
ire(t < t c )/{^e{tc)) « 0. Therefore, < 2 (t)/« 2 (t c )) 
term can be omitted in Eq. (fTO)) . The factor 5/4 in 
front of first T/T* term was introduced for convenience 
purpose. The integration over T\ gives 



S 



£1 
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T 

T* 



cfT_ 
2 T* 
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with qualitative experimental temperature dependencies 
of S in under-doping and near optimal doping. 

In the Eqs. ((9]) and (fTTj) we have used by hand in- 
troduced linear power dependencies on t/t c and T/T*, 
respectively. Further microscopical calculation can es- 
tablish their exact law. However, the goal of the present 
paper was to show that at small these parameters domi- 
nant becomes the single particle boson contribution into 
low temperature heat conductivity and entropy, which is 
probably consistent with experiment. 

As summary we have shown that observed normal state 
low-temperature dependencies of the heat conductivity 
k ~ T 3 - 6 and k ~ T and entropy S - T n with n > 1 
for under-doped and S ~ T for optimal doped cuprate 
superconductors could be the result of functions c ~ T 4 
and c ~ T of the specific heat c. The temperature de- 
pendencies of c were obtained for mixture of IGBQ of 2D 
Coulomb Bose gas and gas of 2D LFLT quasi-particles. 
The obtained from this mixture Lorentz ratio values of 
WFL were close to observed non- Fermi liquid ones. The 
vanishing of free quasi-particles of IGBQ in the transition 
to plasmons has microscopically clarified the decoupling 
of charge carriers with heat and thus, the downturn of 
heat conductivity k. 
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